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Abstract 



Cross sections for the interactions of ultrahigh-energy neutrinos with nucleons are 
evaluated in light of new information about nucleon structure functions. For 10 20 -eV 
neutrinos, the cross section is about 2.4 times previous estimates. We also review the 
cross sections for neutrino interactions with atomic electrons. Some consequences 
for interaction rates in the Earth and for event rates from generic astrophysical 
sources in large-scale detectors are noted. 
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1 Introduction 



Neutrino telescopes hold great promise for probing the deepest reaches of 
stars and galaxies [1-4]. As highly stable neutral particles, neutrinos arrive 
at a detector on a direct line from their source, undeflected by intervening 
magnetic fields. Whereas high-energy photons are completely absorbed by a 
few hundred grams/cm 2 of material, the interaction length of a 1-TeV neutrino 
is about 250 kilotonnes/cm 2 , which corresponds to a column of water 2.5 
million kilometers deep. The feebleness of neutrino interactions means that 
neutrinos can bring us astrophysical information that other radiation cannot, 
but it also means that vast detectors are required to receive this information. 

Encouragement to contemplate neutrino telescopes with effective volumes as 
large as 1 km 3 comes from the observation of neutrinos correlated with su- 
pernova SN1987A [5] and from the detection of solar neutrinos not only by 
radiochemical methods [6-8] but also by observing the direction of recoil elec- 
trons from neutrino interactions [9] . At the same time, detection of neutrinos 
produced by cosmic-ray interactions in Earth's atmosphere [10,11] has be- 
come commonplace in underground detectors [12] and has emerged as a tool 
for investigating neutrino oscillations [13-15]. 

A principal scientific goal of large-scale neutrino telescopes is the detection 
of ultrahigh-energy (UHE: ^ 10 12 eV) cosmic neutrinos produced outside the 
atmosphere: neutrinos produced by galactic cosmic rays interacting with in- 
terstellar gas, and extragalactic neutrinos [16,17]. Extragalactic sources range 
from the conventional — the diffuse (~ 10 18 eV) neutrino flux produced by 
interactions over cosmological time of extragalactic cosmic rays with the mi- 
crowave background radiation [18] — to the highly speculative — such as the 
diffuse flux associated with the decay of cosmic strings [19,20] and other topo- 
logical defects [21] in the relatively late Universe . 

Active galactic nuclei (AGNs) have long been considered as prodigious par- 
ticle accelerators [22] and beam dumps [23], for they are the most powerful 
radiation sources known in the Universe, with typical luminosities in the range 
10 42 to 10 48 erg/s. These cosmic accelerators are presumably powered by the 
gravitational energy of matter spiraling in to a supermassive (~ 1O 8 M ) black 
hole. Cosmic rays generated within an AGN may interact with matter or ra- 
diation in the AGN accretion disk, or with UV photons in the associated jets, 
to produce pions whose decay products include photons and neutrinos. The 
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dominant mechanisms for photon and neutrino production are 



v ip/i) 



7T° + anything 



(1) 




and 



P ip/l) — ^ 7r ± + anything 
U ez/ e z/„ . 



(2) 



If 7r + , 7r~, and 7r are produced in equal numbers, the relative populations 
of the neutral particles will be 27 : 2z/ M : 2z/ M : lu e : lu e . Taken together, 
neutrino emission from ordinary AGNs may provide the dominant isotropic 
flux at energies above about 10 4 GeV. 

The recent detection of energetic photons (E 1 > 100 MeV) from some 40 
AGNs in the Energetic Gamma- Ray Experiment Telescope (EGRET) full-sky 
survey [24] may signal the existence of individual point-sources of neutrinos. 
EGRET, a multilevel thin-plate spark chamber device aboard the Compton 
Gamma-Ray Observatory, has also detected more than a dozen extragalac- 
tic sources at photon energies above 1 GeV. The EGRET sources have the 
characteristics of blazars, AGNs that have associated jets closely aligned with 
the observer's line of sight. The closest EGRET source is Markarian 421, a 
BL Lacertae object at redshift z = 0.031. In 1992, Mrk 421 was detected in 
air showers as a source of TeV photons by the ground-based Whipple Ob- 
servatory, an optical reflector with a 10-meter aperture viewed by more than 
100 small phototubes [25]. In 1995, the Whipple Observatory Gamma- Ray 
Collaboration detected a second TeV photon source, Mrk 501, at z = 0.034 
[26]. If the TeV photons are products of tt° decay, then these sources should 
also be copious neutrino emitters. If instead the TeV photons are produced 
by inverse Compton scattering of energetic electrons off ultraviolet photons, 
no UHE neutrinos will be created. The ability to observe UHE neutrinos from 
TeV photon sources would be an important new AGN diagnostic. 

Ultrahigh-energy neutrinos can be detected by observing long-range muons 
produced in charged- current neutrino-nucleon interactions. To reduce the back- 
ground from muons produced in the atmosphere, it is advantageous to site a 
neutrino telescope at a depth of several kilometers (water equivalent) or to 
observe upward-going muons. High neutrino energy brings a number of ad- 
vantages. First, the charged-current cross section increases, as a oc E u for 
E v £ 10 12 eV, then as a oc E° U A for E v £ 10 15 eV. Second, the background 
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of atmospheric neutrinos falls away compared to the flux from extragalac- 
tic sources, approximately as E~ 1S . Cosmic neutrinos reflect the cosmic-ray 
spectrum near the source (dN/dE oc E~ 2 ), whereas the atmospheric neutrino 
spectrum (oc E~ 3 - 6 above 100 GeV) is about one power of the energy steeper 
than the cosmic-ray spectrum at the Earth (oc E~ 2 - 7 ), which is steeper than 
the source spectrum [27]. The signal of interest for neutrino astronomy should 
emerge from the atmospheric-neutrino background at E v ~ 1-10 TeV. Third, 
the muon range grows with energy, increasing as E^ for E^^l TeV, then in- 
creasing roughly as \ogE^ at higher energies. For upward-going muons, the 
effective volume of a neutrino telescope is thus equal to the instrumented area 
times the muon range. 

Estimates of the fluxes of UHE neutrinos from AGNs and other astrophysi- 
cal sources suggest that a surface area exceeding 0.1 km 2 is required [28]. If 
the muons are detected by observing the Cerenkov light they produce when 
traversing a transparent medium of water or ice, huge target volumes are con- 
ceivable [29]. Four instruments specifically designed for high-energy neutrino 
detection are currently under construction: DUMAND [30,31], at a depth of 
4760 m in the ocean 30 km off the island of Hawaii; the Baikal Neutrino 
Telescope [32], at a depth of 1 km in Lake Baikal in Siberia; NESTOR [33], 
3500 m deep in the Mediterranean near Pylos, Greece; and AMANDA [34,31], 
in deep polar ice at the South Pole. All these detectors aim for effective areas 
of about 0.02 km 2 and an angular resolution for TeV muons of approximately 
1°. These detectors represent a giant step in instrumented volume from their 
underground predecessors. To reach an effective volume of 1 km 3 will require 
efficiencies of scale for the water-Cerenkov technique or new means of de- 
tection. Radio detection is under active study [35]. Acoustic detection may 
become viable in the future [36]. 

At low neutrino energies (E v <C M^/2M, where M w is the intermediate- 
boson mass and M is the nucleon mass), differential and total cross sections 
for the reaction vN — > /j, + anything are proportional to the neutrino energy. 
Above E v ps 10 12 eV, the gauge-boson propagator restricts the momentum 
transfer Q 2 to values near and damps the cross section. At ultrahigh 
energies, the W propagator limits the effective interval in the fractional parton 
momentum x to the region around M^/2ME U . 

Since the UHE vN cross sections were studied in detail nearly a decade ago 
[37-40], our knowledge of parton distributions has developed significantly. In 
place of parton distributions that were essentially based on a single data set 
[41], we now have at our disposal a number of sets of parton distributions 
derived from global fits to a rich universe of experimental information. At small 
values of x, parton distributions have been shaped by measurements made 
possible for the first time by the electron-proton collider HERA at DESY. 
The discovery of the top quark [42,43] with a mass m t ps 175 GeV/c 2 reduces 
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the contribution of the 6-quark sea to the neutrino-nucleon total cross section. 
This new information provides the incentive to reexamine the cross sections 
for UHE vN interactions [44]. 

In §2, we review what is known about the structure of the nucleon and explain 
how we treat the extrapolation to small values of x that is crucial at the 
highest energies. Then in §3 we present in turn our calculations of the charged- 
current and neutral-current cross sections, and explore the variations due to 
different sets of parton distributions. Although vN interactions provide the 
dominant signal and account for most of the attenuation of neutrino beams 
in the Earth at high energies, the W~ resonance in the v e e channel has a very 
strong effect for neutrino energies around 6.3 PeV. Accordingly, we review the 
interactions of neutrinos with electron targets in §4. Section 5 is devoted to a 
study of the attenuation of neutrinos in the Earth. We improve our treatment 
of this important effect by using a detailed model of the Earth's interior. 
We make some remarks about neutrino interactions in the atmosphere in §6, 
and comment in §7 on the possibility of observing the shadows of the Moon 
and Sun. In §8 we estimate the event rates from atmospheric neutrinos and 
from a variety of astrophysical sources in detectors with effective volumes of 
0.1-1 km 3 . A final assessment concludes the paper. 

We find that current knowledge of the proton's parton distributions allows 
us to calculate the uN cross sections with confidence up to neutrino ener- 
gies of about 10 16 eV. The new cross sections are noticeably larger than those 
calculated a decade ago for energies above about 10 15 eV. At 10 20 eV, our 
nominal cross sections are about 2.4 times as large as those calculated using 
the EHLQ parton distributions [37,38]. At energies exceeding 10 16 eV, our ig- 
norance of proton structure at small values of x is reflected in a spread of the 
cross sections calculated using various modern parton distributions. The re- 
sulting uncertainty reaches a factor of 2 ±1 at 10 20 eV. The larger cross sections 
imply enhanced rates for downward-going muons produced in charged-current 
interactions. At the energies of interest for the observation of extraterrestrial 
neutrino sources, upward-going muon rates are little changed, because the 
increased reaction rate is compensated by increased attenuation of neutri- 
nos traversing the Earth to reach the detector. We find that a detector with 
an effective area of 0.1 km and a muon energy threshold in the range of 1 
to 10 TeV should readily observe the diffuse flux of neutrinos expected from 
AGNs above the background of atmospheric neutrinos. The detection of cos- 
mic neutrinos from the interaction of cosmic-ray protons with the microwave 
background appears a remote possibility, even for a 1-km 3 detector. 
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2 New Information about Nucleon Structure 



To compute the cross sections for neutrino-nucleon interactions at high en- 
ergies, we require both a knowledge of the elementary matrix elements and 
also a detailed description of the quark structure of the nucleon. We have 
the first, thanks to extensive experimental validation of the SU(2)l <8> U(l)y 
electroweak theory and refinement of the parameters that appear in the ele- 
mentary neutrino-quark scattering. For the second, we rely on parton distri- 
bution functions extracted from studies of lepton-hadron scattering and of the 
productions of jets, intermediate bosons, dileptons, and photons in hadron- 
hadron collisions. Systematic global fits to experimental data have greatly 
extended our knowledge of parton distribution functions and made modern 
parametrizations increasingly robust. 

Although many experiments have nourished the steady improvement of the 
parton distributions, recent results from the ep collider HERA [45-49] are par- 
ticularly informative for the application at hand. Measurements by the ZEUS 
and HI collaborations mark the first experimental studies of very small parton 
momentum fractions x at momentum transfers Q 2 securely in the deeply in- 
elastic regime. The HERA experiments have begun to map the structure func- 
tion F 2 (x,Q 2 ) in the interval KT 4 &x £ 1(T 2 , with 8.5 GeV 2 £Q 2 £15 GeV 2 . 
For i^2x 10~ 2 , F 2 has been measured over a significant range in Q 2 . 

For most hard-scattering applications in particle physics, it is straightforward 
to begin with parametrizations of parton distribution functions tied to data 
at modest values of Q 2 and evolve them to the desired high scale using the 
Altarelli-Parisi equations [50]. The special challenge of UHE neutrino-nucleon 
scattering is that the PU-boson propagator emphasizes smaller and smaller 
values of x as the neutrino energy E u increases. In the UHE domain, the most 
important contributions to the vN cross section come from x ~ M^/2ME U . 
Up to E v m 10 5 GeV, the parton distributions are sampled only at values of x 
where they have been constrained by experiment. At still higher energies, we 
require parton distributions at such small values of x that direct experimental 
constraints are not available, not even at low values of Q 2 . 

The theoretical uncertainties that enter the evaluation of the UHE neutrino- 
nucleon cross section arise from the low-Q 2 parametrization, the evolution 
of the parton distribution functions to large values of Q 2 ~ Myy, and the 
extrapolation to small values of x. The greatest uncertainty is due to the 
small-x extrapolation. 

Because experiments are limited to values of x ^ 10~ 4 , fits to structure func- 
tions have to be based on plausible but poorly constrained extrapolations to 
x — 0. The parton distributions are traditionally obtained by assuming com- 
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pact forms at Q 2 = Ql = a few GeV 2 : 



xq v (x, Qq) 
xq s (x,Ql) 
xG(x,Ql) 



A v x^(l 
A s x~\l 
A g x-\1 



xrMVx') , 

x) V9 f 9 (Vx~) , 



(3) 



where q v is a valence-quark distribution, q s is a sea-quark distribution, and 
G{x) is the gluon distribution. The functions fi(\/x) are polynomials in -Jx 
that satisfy /i(0) = 1. Sum rules provide broad constraints on the parame- 
ters. For example, the requirement that the momentum integral of the gluon 
distribution be finite means that xG(x, Q 2 ) must be less singular than x" 1 at 
x = 0. The parameters are determined from fits to experimental data and the 
resulting forms are evolved to higher values of Q 2 using the next-to-leading 
order Altarelli-Parisi equations. We employ in this work the latest (CTEQ3) 
of the parton distributions determined by the CTEQ collaboration [51] and 
several sets (MRS A', G, D_, and D_') from the family of parton distributions 
produced by Martin, Roberts, and Stirling [52-55] . Both the CTEQ and MRS 
parametrizations result from global fits to vast data sets and obey sum-rule 
constraints. 

Currently there are two theoretical approaches, both based on perturbative 
QCD, to understanding the (^-evolution of small-rc parton distributions. The 
traditional approach, followed in the CTEQ3 [51] and in the MRS A' [52] and 
G [53] distributions, is to determine parton densities for Q 2 > Ql by solving 
the next-to-leading-order Altarelli-Parisi evolution equations numerically. The 
second approach to small-a: evolution is to solve the Balitskh-Fadin-Kuraev- 
Lipatov (BFKL) equation, which is effectively a leading a s ln(l/:r) resumma- 
tion of soft gluon emissions [56]. In practical terms, the small- a; behavior is an 
input at some scale Ql in the traditional approach, and a dynamically gen- 
erated output in the BFKL scheme. The BFKL approach predicts a singular 
behavior in x and a rapid Q 2 -variation, 



Applying the Altarelli-Parisi equations to singular input distributions oc x 2 
leads to 



a less rapid growth with Q 2 . 

The Altarelli-Parisi approach is applicable in the not-so-small-x and large- 
Q 2 region, while the BFKL solution applies to the small-a; and moderate-Q 2 




(4) 



xq s (x, Q 2 ) ~ In (Q 2 )x 



-0.5 



(5) 
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region. BFKL evolution eventually breaks down at large Q 2 , because of the 
rapid growth exhibited in (4). In the case of ultrahigh-energy neutrino-nucleon 
interactions, the region of interest is small- a: and large-Q 2 , which requires 
a resummation of both \nl/x and \nQ 2 /Ql contributions. Although some 
progress has been made in developing a "unified" evolution equation [57,58], 
the full solution and global fits to data are far from being achieved. 

The standard Altarelli-Parisi evolution of the parton distribution functions is 
applicable for the calculation of the total neutrino cross section up to E v 10 5 
GeV, so it is a reasonable starting point for calculating the cross section for 
higher energy neutrinos. Consequently, the calculation of the total neutrino- 
nucleon cross section presented here relies on the CTEQ3 and MRS A' parton 
distributions obtained using next-to-leading-order (NLO) evolution equations. 
The CTEQ3 distribution functions, depending on the order of the evolution 
and the factorization scheme, use A ~ 0.28 — 0.35, while A = 0.17 for MRS A'. 
The CTEQ3 distributions are particularly convenient as a benchmark because 
the numerical evolution is provided for x — > 0, including the region in which 
the Altarelli-Parisi equations may not be reliable. (The MRS A' distributions 
are available for x > 10~ 5 .) We use the CTEQ3 parton distributions, with 
NLO evolution from Q = 1.6 GeV, as our canonical set. We choose the deep- 
inelastic scattering factorization scheme (DIS) parametrization of the parton 
distribution functions, for which the exponent A = 0.332. Results calculated 
with this set of parton distributions are labeled as CTEQ-DIS in the discussion 
below. 

To estimate the uncertainty in the small-x parton distributions evaluated 
at Q 2 ~ M$y, we consider alternative treatments of the small-x behavior. 
To explore a less singular alternative, we extrapolate to x = using the 
double-logarithmic-approximation (DLA) [59] , an approximate solution to the 
Altarelli-Parisi equations for not-too-singular input distributions. The form of 
the sea-quark distribution is [59,60] 



Here, p = (8N/b )\n(l/x), f(Q) = lnln(Q 2 /A 2 ), A = 3 is the number of 
colors, and 6o — (HA— 2n/)/3 forn/ flavors. This form was used in Ref. [37,38] 
to extrapolate the EHLQ parton distributions [41] below x mm = 10~ 4 . The 
EHLQ distributions xq s (x, Qq) are finite as x — > 0, i.e., correspond to A = 
in (4) [61]. To estimate a lower limit on the UHE vN cross section, we use the 
DLA form of Eq. (6) for x < x mm = 10~ 4 as an extrapolation of the leading- 
order parametrization of the CTEQ3 parton distribution functions (labeled 
CTEQ-DLA). Here, following the procedure of Ref. [37], we choose C(Q 2 ) 



to match x min q s (x mh \Q 2 ) from CTEQ-LO. The five-flavor value of A QCD is 
Agg D = 132MeV[62]. 




(6) 
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X 

Fig. 1. Comparison of the light-quark sea at Q 2 = M^y for various parton distribu- 
tions. Of the MRS distributions, D_ (A') is the most (least) singular. 

A more singular form of the parton distributions at small-x, motivated by 
BFKL dynamics, appears in the MRS D_ set [54]. In the limit of very small 
x, the behavior of the MRS D_ sea is xq s (x,Ql) = C(Ql)x~ ' 5 . These dis- 
tributions appear to slightly overestimate the low-Q 2 HERA data [48,49] in 
the interval 1CT 4 ^x ^ 1CT 2 . Thus, for large E u , they can provide a reasonable 
upper limit on the cross section [63]. 

To illustrate the range of parton distributions that these choices represent, 
we plot the light-quark sea distribution x(u + d)/2 versus x for Q 2 = 
in Figure 1. For a;^10~ 4 , the MRS and CTEQ distributions all are in close 
agreement. This consonance allows us to make confident predictions of the uN 
cross sections for neutrino energies up to about 10 6 GeV. The spread in the 
parton distributions at smaller values of x reflects the uncertain extrapolation 
toward x — 0. 



3 Neutrino— Nucleon Interactions 

It is straightforward to calculate the inclusive cross section for the reaction 

v^N -> fT + anything, (7) 

n + p 

where N = — - — is an isoscalar nucleon, in the renormalization group- 
improved parton model. The differential cross section is written in terms of 
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the Bjorken scaling variables x = Q 2 /2Mv and y = v/E u as 

72, o^<9 7. r 7-1 / 7. T \ 2 



d 2 a 2G 2 F ME U 



M 2 W 



dxdy 



TT \Q 2 + M 2 V/ 



xq(x,Q 2 )+xq(x,Q 2 )(l-y) 2 , (8) 



where — Q 2 is the invariant momentum transfer between the incident neutrino 
and outgoing muon, v = E v —E^ is the energy loss in the lab (target) frame, M 
and Mw are the nucleon and intermediate-boson masses, and Gf = 1.16632 x 
1CT 5 GeV~ 2 is the Fermi constant. The quark distribution functions are 

/ u v (x,Q 2 ) + d v (x,Q 2 ) u s (x,Q 2 ) +d s (x,Q 2 ) 

Q ) = 2 + 2 

+s s (x,Q 2 ) + b s (x,Q 2 ) (9) 

q(x, Q 2 ) = u ^Q 2 ) + dsM2) + c s (x, Q 2 ) + t s (x, Q 2 ), 

where the subscripts v and s label valence and sea contributions, and u, d, c, 
s, t, b denote the distributions for various quark flavors in a proton. At the 
energies of interest for neutrino astronomy, perturbative QCD corrections to 
the cross section formula (8) are insignificant, so we omit them. In particular, 
in the DIS factorization scheme (the CTEQ-DIS parton distributions), the 
terms proportional to a s [65] in the NLO cross section contribute only a few 
percent. 

Because of the great mass of the top quark, ti pairs are a negligible component 
of the nucleon over the Q 2 -range relevant to neutrino-nucleon scattering. Con- 
sequently we drop the contribution of the top sea. At the energies of interest 
here, it is a sound kinematical simplification to treat charm and bottom quarks 
as massless. However, the threshold suppression of the b — > t transition must 
be taken into account. We adopt the standard "slow-rescaling" prescription 
[64], with m t = 175 GeV/c 2 . Numerical integrations were carried out using the 
adaptive Monte Carlo routine VEGAS [66], and Gaussian techniques. 

We show in Figure 2 the contributions of valence quarks and of the different 
quark flavors in the sea to the vN charged-current total cross section, accord- 
ing to the CTEQ3 parton distributions. As expected, the valence contribution 
dominates at low energies. There, in the parton-model idealization that quark 
distributions are independent of Q 2 , differential and total cross sections are 
proportional to the neutrino energy. Up to energies E u ~ 10 11 eV, the fa- 
miliar manifestation of the QCD evolution of the parton distributions is to 
degrade the valence component, and so to decrease the total cross section. At 
still higher energies, the gauge-boson propagator restricts Q 2 = 2ME u xy to 
values near M^, and so limits the effective interval in x to the region around 
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Fig. 2. Components of the vN charged-current cross section as functions of the 
neutrino energy for the CTEQ3 distributions. 



M^ v /2ME u (y) . Figure 3 shows the contributions to the cross section from dif- 
ferent regions of x. At modest values of Q 2 , the effect of this PF-propagator 
damping is to further diminish the cross section below the point-coupling, 
parton-model approximation. Above about 10 16 eV, the valence contribution 
is even smaller than the contribution of the bb sea. 

A second effect of QCD evolution is to increase the population of heavy quarks 
(s, c, b) within the proton, and to increase the importance of the light-quark 
sea at small values of x. Andreev, Berezinsky, and Smirnov [67] have pointed 
out that the effect of this growth in the density of the parton sea is to enhance 
the cross section at high energies. This effect is apparent in Figure 4. There we 
compare the CTEQ3 cross section with the 1986 cross section [37] based on the 
EHLQ structure functions and with the case of no evolution. We see that the 
EHLQ-based cross section is enhanced by fully an order of magnitude at high 
energies by the evolution of the sea. At low energies, the decrease in the cross 
section brought about by the degradation of the valence distribution is appar- 
ent in the comparison of the EHLQ curves with and without evolution. We 
also show in Figure 4 the CTEQ3 prediction for the vN charged-current cross 
section. At the highest energies, where the contributions of valence quarks are 
unimportant, the neutrino and antineutrino cross sections are identical. 
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Our new evaluation of the vN cross section differs from the earlier calculations 
at both low and high energies. At both extremes, the difference is owed to 
changes in our understanding of parton distribution functions. The EHLQ 
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Fig. 3. Integral cross section (l/a) Jq dxda/dx for the charged-current reaction 
u^N -> fj,~ + anything at E v = 10 5 ,10 7 , and 10 9 GeV. As the neutrino energy 
increases, the dominant contributions come from smaller values of x. 
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cording to the CTEQ3 parton distributions. The EHLQ-DLA prediction [37] for 
the uN cross section is also shown, together with the vN cross section based on the 
unevolved EHLQ structure functions with Q 2 fixed at Qq = 5 GeV 2 . 
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parton distributions, on which the earlier calculations were based, were based 
on the CERN-Dortmund-Heidelberg-Saclay measurements of neutrino-nucleon 
structure functions [68]. We now know that the normalization of the CDHS 
structure functions was about 15% low [69]. The change in normalization 
directly affects the cross sections at low energies. At higher energies, which 
are sensitive to small values of x, the shape of the parton distribution as 
x — > is decisive. At low Q 2 , the EHLQ distributions xq s (x) are finite as 
x — »• 0, whereas HERA experiments point to singular behavior, parametrized 
in the CTEQ distributions as xq s (x) — > a; -0 - 332 . The density of partons at 
small values of x and modest values of Q 2 is thus greater than was assumed 
in the earlier work. 

We show in Figure 5 the charged-current uN cross section implied by sev- 
eral sets of parton distributions derived from global fits. There is excellent 
agreement among the predictions of the MRS D_, G, and A' distributions and 
the CTEQ3 distributions up to E u « 10 7 GeV. Above that energy, our DLA 
modification of the CTEQ3 distributions gives a lower cross section than the 
full CTEQ3 distributions (CTEQ-DIS), as expected from its less singular be- 
havior as x — > 0. At the highest energy displayed, the most singular (MRS 
D_) distribution predicts a significantly higher cross section than the others. 
Above about 10 6 GeV, the EHLQ-DLA distributions yield noticeably smaller 
cross sections than the modern distributions. All the MRS and CTEQ curves 
are in reasonable agreement with the HERA measurement [70] of the charged- 
current cross section at an equivalent neutrino energy of 46.7 TeV [71]. The 
parton distributions inferred by Frichter, et al. from HERA data [44] yield 
cross sections that stand apart from those derived from global fits. 

Two other groups recently have evaluated the neutrino-nucleon charged-current 
cross sections at high energies. Parente and Zas [72] used the MRS G distribu- 
tions [53] to compute a cc (vN) for neutrino energies in the range 200 GeV < 
E u < 10 7 GeV, in which no special treatment of the x — > behavior of the 
parton distributions is required. The results presented in their Figure 2 agree 
with the corresponding curve in Figure 5 above. Butkevich, et al. [73] have 
evaluated a cc (vN) and a cc (vN) for 10 2 GeV < E v < 10 6 GeV using the 
MRS A distributions [52] , an early version of the Gliick-Reya-Vogt distribu- 
tions [74], and the Morfin-Tung ancestor [75] of the CTEQ3 distributions we 
use. The results presented in their Figure 1 agree with those in our Figure 5. 
Butkevich, et al. have also explored two extrapolations of the MRS A parton 
distributions to very small values of x. The values of acciyN) presented in 
their Table 1 for 10 7 GeV < E u < 10 12 GeV are close to those we give in 
Table 1 below for our nominal set, the CTEQ3 distributions, and agree well 
with our calculations using the MRS A' distributions. 

The differential cross section (1 / E u )da / dy for neutrino-nucleon scattering is 
shown in Figure 6. The peaking of the cross section near y — 0, which becomes 
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increasingly prominent with increasing neutrino energy, is a direct consequence 
of the cutoff in Q 2 enforced by the W propagator. However, because of the 
growth of the quark distributions as small values of x for large Q 2 , the cross 
section is nonnegligible at finite values of y. Accordingly, the mean inelasticity 
(y) does not decrease rapidly as the energy increases. This parameter is shown 
for both neutrinos and antineutrinos in Figure 7. 

A parallel calculation leads to the neutral-current cross section. In this case 
the differential cross section for the reaction u^N — > v + anything is given by 



d 2 a G 2 F ME U 



M 2 



dxdy 



2tt \Q 2 + M 2 Z/ 



xq (x,Q 2 )+xq°(x,Q 2 )(l-y) 2 ], (10) 



where M z is the mass of the neutral intermediate boson. The quantities in- 
volving parton distribution functions are 



q°(x,Q 2 ) 



u v (x, Q 2 ) + d v (x, Q 2 ) u s (x, Q 2 ) + d s (x, Q 2 ) 



+ 



(Ll + L 2 ) 



+ 



u s (x, Q 2 ) + d 8 (x, Q' 



[s s (x, Q 2 ) + b s (x, Q 2 )}(L 2 d + R 2 ) + [c s (x, Q 2 ) + t s (x, Q 2 )}(L 2 U + R 2 U ) 
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Fig. 5. The charged-current cross section for v„ interactions with an 
isoscalar nucleon. The parametrization of Frichter, et al. [44] is shown for 
5 x 10 4 GeV < E v < 5 x 10 7 GeV. The data point is an average of measurements 
by the ZEUS and HI Collaborations at HERA [70]. 
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+ 



u v (x, Q 2 ) + d v (x, Q 2 ) u s (x,Q 2 ) + d s (x,Q 



+ 



u s (x,Q 2 ) + d s (x,Q 2 ) 



(Ll + L 2 ) + 



(12) 



[s s (x, Q 2 ) + b s (x, Q 2 )}(L 2 d + R 2 ) + [c s (x, Q 2 ) + t s (x, Q 2 )}(L 2 U + R 2 U ), 





10 100 1000 10 4 10 5 10 6 10 7 10 8 10 9 10 10 10 n 10 12 

E y [GeV] 

Fig. 7. Energy dependence of the inelasticity parameter y for charged-current (solid 
lines) and neutral-current (dashed lines) interactions as a function of the incident 
neutrino energy. 
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where the chiral couplings are 

i-r„, r„ = _u 

(13) 



1 — \xyf Ld — —1 + \xw 



R u — —^xw Rd — \xw 



and xw = sin 2 6w is the weak mixing parameter. For numerical calculations 
we have chosen x w = 0.226 [76]. Again the top-quark sea is negligible. 

Cross sections for neutral-current scattering of neutrinos and antineutrinos 
from isoscalar nucleons are shown as the dashed curves in Figures 8 and 9, 
respectively. There we also show the charged-current cross sections (as thin 
solid curves) and the sum of charged-current and neutral-current cross sections 
(as thick solid curves). 

Numerical values of the cross sections and inelasticity parameters, which char- 
acterize the angular distribution of outgoing leptons, are indispensable for sim- 
ulating the degradation of the neutrino flux passing through the Earth, and 
for calculating event rates in proposed detectors. We have gathered in Tables 
1 and 2 the charged-current and neutral-current cross sections and values of 
(y), for vN and uN collisions, respectively. 

For neutrino energies in the range 10 15 eV < E v < 10 21 eV, good representa- 
tions of the cross sections are given by simple power-law forms: 
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E y [GeV] 

Fig. 8. Cross sections for vN interactions at high energies: dotted line, 
a{vN — > v + anything); thin line, <j(vN — > fi~ + anything); thick line, total 
(charged-current plus neutral-current) cross section. 
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Table 1 

Charged-current and neutral-current cross sections for vN interactions, and the 
corresponding values of the mean inelasticity (y), for the CTEQ-DIS distributions. 



E v [GeV] 


ace [cm 2 


] 


cr NC [cm 2 


] 


(y)cc 


{yhc 


10 1 


0.777 x 10" 


-37 


0.242 x 10" 


-37 


0.483 


0.474 


10 2 


0.697 x 10" 


-36 


0.217 x 10" 


-36 


0.477 


0.470 


10 3 


0.625 x 10" 


-35 


0.199 x 10" 


-35 


0.472 


0.467 


10 4 


0.454 x 10" 


-34 


0.155 x 10- 


-34 


0.426 


0.428 


10 5 


0.196 x 10" 


-33 


0.745 x 10" 


-34 


0.332 


0.341 


10 6 


0.611 x 10- 


-33 


0.252 x 10- 


-33 


0.273 


0.279 


10 7 


0.176 x 10" 


-32 


0.748 x 10- 


-33 


0.250 


0.254 


10 8 


0.478 x 10" 


-32 


0.207 x 10" 


-32 


0.237 


0.239 


10 9 


0.123 x 10- 


-31 


0.540 x 10" 


-32 


0.225 


0.227 


10 10 


0.301 x 10" 


-31 


0.134 x 10" 


-31 


0.216 


0.217 


10 11 


0.706 x 10" 


-31 


0.316 x 10" 


-31 


0.208 


0.210 


10 12 


0.159 x 10" 


-30 


0.715 x 10" 


-31 


0.205 


0.207 



Table 2 

Charged-current and neutral-current cross sections for uN interactions, and the 
corresponding values of the mean inelasticity (y), for the CTEQ-DIS distributions. 

E v [GeV] ace [cm 2 ] cr NC [cm 2 ] (y) C c (j/)nc 



10 1 


0.368 x 10- 


-37 


0.130 x 10- 


-37 


0.333 


0.350 


10 2 


0.349 x 10' 


-36 


0.122 x 10" 


-36 


0.340 


0.354 


10 3 


0.338 x 10" 


-35 


0.120 x 10" 


-35 


0.354 


0.368 


10 4 


0.292 x 10" 


-34 


0.106 x 10" 


-34 


0.345 


0.358 


10 5 


0.162 x 10" 


-33 


0.631 x 10- 


-34 


0.301 


0.313 


10 6 


0.582 x 10- 


-33 


0.241 x 10" 


-33 


0.266 


0.273 


10 7 


0.174 x 10" 


-32 


0.742 x 10- 


-33 


0.249 


0.253 


10 8 


0.477 x 10- 


-32 


0.207 x 10" 


-32 


0.237 


0.239 


10 9 


0.123 x 10" 


-31 


0.540 x 10" 


-32 


0.225 


0.227 


10 10 


0.301 x 10" 


-31 


0.134 x 10" 


-31 


0.216 


0.217 


10 11 


0.706 x 10- 


-31 


0.316 x 10- 


-31 


0.208 


0.210 


10 12 


0.159 x 10" 


-30 


0.715 x 10" 


-31 


0.205 


0.207 
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Fig. 9. Cross sections for vN interactions at high energies: dotted line, 
a{vN — > v + anything); thin line, a{pN — > /x + + anything); thick line, total 
(charged-current plus neutral-current) cross section. 

/ T? \ 0-402 

<rcc(fW) =2.69 x 10- 36 cm 2 (^^J 



a NC (z/AT) = 1.06 x 1CT 36 cm 2 
<7 CC (pjV)=2.53 x 10" 36 cm 2 
a NC (z/Ar)=0.98 x 1CT 36 cm 2 



E v \°- 408 
1 GeV J 

E v \ 0404 
1 GeV J 

1 GeV/ 



(14) 



0.410 



Before leaving the subject of neutrino- nucleoli collisions, let us note that the 
cross sections for the reactions vN W oy Z + anything are small compared 
with the cross sections for deeply inelastic scattering [77]. 



4 Interaction of UHE Neutrinos with Electrons 



Because of the electron's small mass, neutrino-electron interactions can gen- 
erally be neglected with respect to neutrino-nucleon interactions [78]. There is 
one exceptional case: resonant formation of the intermediate boson W~ in v e e 
interactions at 6.3 PeV [79]. The resonant cross section is larger than the vN 
cross section at any energy up to 10 21 eV. Accordingly, it is important to have 
the neutrino-electron cross sections in mind when assessing the capabilities of 
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neutrino telescopes. 



Defining as usual the laboratory energy of the incoming neutrino as E v and 
the laboratory energy of the recoiling charged lepton as E' = yE u , we may 
write the differential cross sections for neutrino-electron scattering as [80] 



daiy^e — > v^e) _ G F mE u 



dy 2tt (i + 2mE v y/M 2 ) 

daip^e — > P M e) _ G F mE u 1 



Rl{l-yf + Ll 
R 2 e + L 2 e (l-y) 2 



dy 2tt (i + 2mE u y/M Z/ 

doiv^e — >• //i/ e ) _ G 2 F mE v 4[1 - - m 2 )/2mE i/ ] 2 



2vr (l + 2m^(l-y)/M^) 



2 \2 ' 



(15) 
(16) 
(17) 



da(y e e — >■ z/ e e) G 2 F mE v 



dy 



2tt 



i£(l-v) 5 



(1 + 2mE v y/M 2 z f 



+ 



+ 



1 + 2mE v y/M 2 1 + 2m£„(l - 



(18) 



da(u e e — > z/ e e) G 2 F mE u 



dy 



2n 



+ 



;i + 2mE v y/M 2 ) 2 
2 



+ 



1 + 2mE v y/M 2 z 1 - 2mE u /M$ v + iT w /M w 



da(v e e -> z^/i) _ G 2 F mE u 4(1 - y) 2 [l - (^ 2 - m 2 )/2mKf 



2tt (l-2mE I/ /M^ + ryM 2 , ' 



(19) 



(20) 



and 



da{y e e — > hadrons) da(u e e — > z^/i) T(W^ — > hadrons) 



r(w - //p M ) 



(21) 



where m = 0.51099908 MeV/c 2 is the electron mass and // = 105.658389 MeV/c 2 
is the muon mass [81]. The chiral couplings of the Z° to the electron are 
L e = 2sm 2 9 w - 1 and R e = 2sin 2 ^, with sin 2 6 W = 0.226 [76]. To evalu- 
ate the cross sections, we use M w = 80.22 GeV/c 2 , M z = 90.188 GeV/c 2 , and 
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Fig. 10. Cross sections for neutrino interactions on electron targets. At low energies, 
from largest to smallest cross section, the processes are (i) v e e — ► hadrons, (ii) 
u^e -> five, (iii) v e e -»■ f e e, (iv) P e e -> z^//, (v) v e e -»■ P e e, (vi) f M e -»■ i/^e, (vii) 



are plotted in Figure 10. Only in the neighborhood of the intermediate-boson 
resonance are any of the neutrino-electron processes competitive with the 
neutrino-nucleon cross sections. The cross sections at the resonance peak, 
E r u cs = Myy/Zrn, are collected in Table 3, together with the cross sections 
for neutrino-nucleon scattering. 

We shall consider the effects of the W~ resonance region, (M w — 2T w ) 2 /2m = 
5.7 PeV & E v &(M W + 2T w ) 2 /2m = 7.0 PeV, on the attenuation of cosmic v e 
in the Earth, through the reaction v e e — > W~ — > anything, in §5. In §8.3 
we project the rate of downward-going v e e — > W~ — > u^fi events for various 
models of the diffuse neutrino flux from active galactic nuclei. 



Y w = 2.08 GeV. The integrated cross sections 




(22) 
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Table 3 

Integrated cross sections for neutrino-electron and neutrino-nucleon scattering at 
El cs = M^/2m = 6. 3 x 10 6 GeV. 

Reaction a [cm 2 ] 





-»■ v^e 


5. 


.86 


X 


10 


-OO 


v^e 


-> P M e 


5 


.16 


X 


10 


-36 


Vp,e 




5. 


.42 


X 


10 


-35 


v e e 


-»■ v e e 


3 


.10 


X 


10 


-35 


v e e 


-»■ z/ e e 


5 


.38 


X 


10 


-32 


v e e 




5 


.38 


X 


10 


-32 


v e e 


— > P T T 


5 


,38 


X 


10 


-32 


v e e -» 


hadrons 


3 


.41 


X 


10 


-31 


v e e -»■ 


anything 


5. 


.02 


X 


10 


-31 



v^N — ► n~ + anything 1.43 x 10 _: 

u^N —> + anything 6.04 x 10~' 

u^N — > /i+ + anything 1.41 x 10 _: 

z^iV — > P M + anything 5.98 x 10 _: 

5 The Earth is Opaque to UHE Neutrinos 



The rise of the charged-current and neutral-current cross sections with energy 
is mirrored in the decrease of the (water-equivalent) interaction length, 

£mt = ^bvA ' (23) 



where N\ = 6.022 x 10 23 mol -1 = 6.022 x 10 23 cm -3 (water equivalent) is 
Avogadro's number. The energy dependence of the interaction lengths for neu- 
trinos on nucleons is shown in Figure 11. We show separately the interaction 
lengths for charged-current and neutral-current reactions, as well as the inter- 
action length corresponding to the total (charged-current plus neutral-current) 
cross section. The same information is shown for antineutrinos on nucleons in 
Figure 12. Above about 10 16 eV, the two sets of interaction lengths coincide. 
These results apply equally to u e N (or u e N) collisions as to u^N (or v^N) 
collisions. 

Over the energy range of interest for neutrino astronomy, the interactions of v e , 
Vp, and P M with electrons in the Earth can generally be neglected in compari- 
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E v [GeV] 

Fig. 11. Interaction lengths for neutrino interactions on nucleon targets: dotted line, 
charged-current interaction length; dashed line, neutral-current interaction length; 
solid line, total interaction length, all computed with the CTEQ-DIS parton distri- 
butions. The dot-dashed curve shows the charged-current interaction length based 
on the EHLQ structure functions with Q 2 fixed at Qq = 5 GeV 2 , as in Figure 4. 

son to interactions with nucleons. The case of u e e interactions is exceptional, 
because of the intermediate-boson resonance formed in the neighborhood of 
El es = M^r /2m « 6.3 x 10 15 eV. The resonant reactions v e e -> W~ -> 
and v e e — * W~ — > hadrons may offer a detectable signal. At resonance, the re- 
action v e e — > W~ — > anything significantly attenuates a v e beam propagating 
through the Earth. The water-equivalent interaction lengths corresponding to 
the neutrino-electron cross sections computed in §4 are displayed in Figure 13. 
These are evaluated as 



^ e (^)(10/18)7V A ' 



where (10/18)A^a is the number of electrons in a mole of water. 

To good approximation, the Earth may be regarded as a spherically symmetric 
ball with a complex internal structure consisting of a dense inner and outer 
core and a lower mantle of medium density, covered by a transition zone, lid, 
crust, and oceans [82]. A convenient representation of the density profile of 
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E v [GeV] 

Fig. 12. Interaction lengths for antineutrino interactions on nucleon targets: dot- 
ted line, charged-current interaction length; dashed line, neutral-current interaction 
length; solid line, total interaction length, all computed with the CTEQ-DIS parton 
distributions. 



the Earth is given by the Preliminary Earth Model [83], 



p(r) 



13.0885 - 8.8381a; 2 , r < 1221.5 

12.5815 - 1.2638a; - 3.6426a; 2 - 5.5281a; 3 , 1221.5 < r < 3480 
7.9565 - 6.4761a; + 5.5283a; 2 - 3.0807x 3 , 3480 < r < 5701 



5.3197- 1.4836a;, 

11.2494 - 8.0298a;, 

7.1089 - 3.8045a;, 

2.691 + 0.6924x, 

2.9, 

2.6, 

1.02, 



5701 < r < 5771 
5771 < r < 5971 
5971 < r < 6151 
6151 < r < 6346.6 
6346.6 < r < 6356 
6356 < r < 6368 
r < i? ffi , 



(25) 



23 




E v [GeV] 

Fig. 13. Interaction lengths for neutrino interactions on electron targets. At low ener- 
gies, from smallest to largest interaction length, the processes are (i) P e e — > hadrons, 
(ii) f M e -»■ (j,i/ e , (hi) v e e -> i/ e e, (iv) P e e -> P^/x, (v) P e e -»■ P e e, (vi) v^e -> i/ M e, (vii) 
z^e -> z^e. 

where the density is measured in g/cm 3 , the distance r from the center of the 
Earth is measured in km and the scaled radial variable x = r/i? e , with the 
Earth's radius R® = 6371 km. The density of a spherically symmetric Earth 
is plotted in Figure 14. 

The amount of material encountered by an upward-going neutrino in its pas- 
sage through the Earth is shown in Figure 15 as a function of the neutrino 
direction. The influence of the core is clearly visible at angles below about 0.2n. 
A neutrino emerging from the nadir has traversed a column whose depth is 11 
kilotonnes/cm 2 , or 1.1 x 10 10 cmwe. The Earth's diameter exceeds the charged- 
current interaction length of neutrinos with energy greater than 40 TeV. In 
the interval 2 x 10 6 GeV ^ E u ^2 x 10 7 GeV, resonant v e e scattering adds 
dramatically to the attenuation of electron antineutrinos. At resonance, the 
interaction length due to the reaction v e e — > W~ — > anything is 6 tonnes/cm 2 , 
or 6 x 10 6 cmwe, or 60 kmwe. The resonance is effectively extinguished for neu- 
trinos that traverse the Earth. 

We discuss the effect of attenuation on interaction rates of upward-going muon- 
neutrinos in §8.2. 
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Fig. 15. Thickness of the Earth as a function of the angle of incidence of the incoming 
neutrinos. 

6 UHE Neutrino Interactions in the Atmosphere 

The atmosphere is more than a thousand times less dense than the Earth's 
interior, so it makes a negligible contribution to the attenuation of the incident 
neutrino flux. The US Standard Atmosphere (1976) [84] can be reproduced to 
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3% approximation by the following simple parametrization: 

1.225 x 1(T 3 g/cm 3 exp (-/i/9.192 km), h< 10 km, 



Patm(^) = < 



(26) 

1.944 x 10^ 3 g/cm 3 exp (-/i/6.452 km) /i > 10 km. 



For a standard atmosphere, a neutrino normally incident on a surface detector 
passes through a column density of 1033 g/cm 2 = 1033 cmwe, while a neu- 
trino arriving along the horizon passes through a column of about 36 000 cmwe. 
Both amounts of matter are orders of magnitude smaller than the neutrino 
interaction lengths at the energies under study (cf. Figures 11, 12, and 13). 
The atmosphere is thus essentially transparent to neutrinos. 

On the other hand, the amount of material encountered by a neutrino passing 
horizontally through the atmosphere is not small compared with the depth 
available for the production of contained events in a water (or ice) Cerenkov 
detector. Figure 16 shows the column depth traversed by a horizontal neutrino 
as a function of altitude. (The values shown are for the full passage through the 
atmosphere, not just inbound to the point of closest approach to the surface.) 
An air shower detector like the Fly's Eye [85] , which detects light produced by 
nitrogen fluorescence along the path of a high-energy particle traversing the 
atmosphere, could detect neutrino-induced cascades and perhaps identify their 
shower profiles. Indeed, Halzen, et al. [86] have argued that the 3 x 10 20 -eV 
cosmic ray shower observed by Fly's Eye [87], the highest energy cosmic-ray 
event, might have been initiated by a neutrino. Sigl and Lee [88] comment that 
the interpretation of the highest-energy cosmic rays as neutrino interactions in 
the atmosphere becomes likely only if <Tcc(vN) were a few orders of magnitude 
higher than we calculate. 



7 Shadows of the Moon and Sun 



In recent years, cosmic-ray experiments have used the observation of shadow- 
ing of the cosmic-ray flux by the Moon and Sun to demonstrate the angular 
resolution of their detectors [89] . Might it someday be possible to observe the 
shadowing of neutrinos by Earth's satellite and star? 

The Moon has a radius of -RMoon = 1738 km and an average density of 
(pMoon) = 3.37 g/cm 3 . It is approximately uniform in density, except for a 
core at R < 238 km, where pMoon ~ 7.55 g/cm 3 [90]. The column depth along 
the lunar diameter is 1.378 x 10 9 cmwe, which makes the Moon opaque to 
neutrinos with E v ^ 10 6 GeV. 

The matter distribution in the Sun extends to a solar radius of R Q = 6.96 x 
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Fig. 16. Column depth encountered by a horizontal neutrino traversing Earth's 
atmosphere at an altitude h. 

10 5 km. The density distribution is known from the standard solar model 
[91]. Except very near the center, a good description is given by the simple 
parametrization, 

p = 236.93 g/cm 3 exp(-10.098 r/R & ). (27) 

The profile through the solar diameter is 3.27 x 10 12 cmwe, which makes the 
Sun opaque to neutrinos with E v ^ 100 GeV. The column density encountered 
by parallel rays of neutrinos falling on the Sun is shown as a function of 
distance from the center of the Sun's face in Figure 17. Almost the entire face 
of the Sun is opaque to neutrinos with energies above 10 6 GeV. 

Since the Moon and Sun are small in the sky, each with an angular diameter 
of about 1/2°, both large detector volumes and excellent angular resolution 
will be required to see their shadows. 



8 UHE Neutrino Fluxes and Event Rates 

In this section, we calculate event rates for atmospheric neutrinos, cosmic 
neutrinos and neutrinos that originate in active galactic nuclei. We start with 
a brief discussion of theoretical models for UHE neutrino production and their 
predictions for the energy dependence of muon-neutrino and electron-neutrino 
fluxes. We consider representative fluxes in order to assess the feasibility of 
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Fig. 17. Column depth presented by the Sun to parallel streams of neutrinos. 

detection and examine the consequences of our new neutrino-nucleon cross 
sections. We first calculate the event rates for upward-moving muons and 
antimuons produced in the material below the detector, and then consider 
rates for downward-moving and contained events for both muon- and electron- 
neutrino interactions. 



8. 1 Sources of UHE Neutrinos 

In Figures 18 and 19 we display differential neutrino fluxes from a variety 
of sources. Neutrinos produced by cosmic-ray interactions in the atmosphere 
dominate other neutrino sources at energies below 1 TeV. For the detection 
of extraterrestrial neutrinos we focus on neutrino energies above 1 TeV. The 
solid curves shown in Figure 18 represent + fluxes produced by several 
mechanisms, while Figure 19 shows the v e + D e fluxes. 

The "conventional" atmospheric neutrino flux at E v — 1 TeV is derived from 
the decay of charged pions and kaons produced by cosmic ray interactions 
in the atmosphere. The conventional flux calculated by Volkova [92], labeled 
by ATM in the figures, is exhibited as the angle average of the atmospheric 
v \i + (Figure 18) and v e + v e (Figure 19) fluxes. The predicted horizontal 
neutrino spectra are in agreement with the absolute spectra measured in the 
Frejus experiment up to 10 TeV [12,93]. The atmospheric neutrino flux is large 
at E v — 1 TeV, but the spectrum falls rapidly as a function of energy. For 
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Fig. 18. Muon neutrino plus antineutrino fluxes at the Earth's surface: an- 
gle-averaged flux from cosmic-ray interactions in the atmosphere (ATM), and 
isotropic fluxes from active galactic nuclei (AGN-SS, AGN-NMB, and AGN-SP) and 
from cosmic-ray interactions with the microwave background (CR-2 and CR-4). The 
Frejus upper limit [101] for a neutrino flux in excess of the atmospheric neutrino 
flux is indicated at 2.6 TeV. The dotted line shows the vertical flux of atmospheric 
/i + + fi~ calculated in Ref. [94]. 

1 TeV < E u < 10 3 TeV, the angle-averaged atmospheric + flux can be 
approximated by a power law spectrum: 



dN Vu 



dE u 



7.8 x 10 



-ii 



TeV/ 



cm 



2 s" 1 sr" 1 GeV" 1 , 



(28) 



The use of the angle-averaged atmospheric flux, while not necessary, facilitates 
comparison with fluxes from diffuse extraterrestrial sources. 

An additional "prompt" contribution to the atmospheric flux arises from 
charm production and decay. The vertical prompt neutrino flux has recently 
been reexamined using the Lund model for particle production [94], and has 
been shown to be small relative to the conventional atmospheric flux for 
E v < 10 5 GeV. Since atmospheric neutrinos are a significant background only 
for E v ^ 10 TeV, we neglect neutrinos from charm decay in our calculations 
of event rates. 

We also show in Figure 18 the vertical atmospheric muon flux from conven- 
tional and prompt sources [94], indicated by a dotted line. The atmospheric 
muon flux for E^ > 10 7 GeV is dominated by muons from charm decays. The 
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Fig. 19. Electron neutrino plus antineutrino fluxes at the Earth's surface: an- 
gle-averaged flux from cosmic-ray interactions in the atmosphere (ATM), and 
isotropic fluxes from active galactic nuclei (AGN-SS, AGN-NMB, and AGN-SP) 
and from cosmic-ray interactions with the microwave background (CR-2 and CR-4). 



muon spectrum at sea level is approximately parametrized by 



dN, 



dE, 



^ - 1.05 x 10- 10 



-3.7 



E H Y -2 - 

ITWJ Cm 8 



1 sr- 1 GeV" 1 . 



(29) 



Atmospheric muons from charm decay and from conventional sources consti- 
tute a background to the detection of v^N charged-current interactions. By 
deploying a detector at great depths [95], or observing upward-going muons, 
or both, one can reduce the cosmic-ray muons to a manageable background. 



Detectable fluxes of neutrinos may be generated in active galactic nuclei [96]. 
The observation [97] that the diffuse neutrino flux from unresolved AGNs 
might be observable with the proposed neutrino telescopes has stimulated a 
number of calculations of the diffuse UHE neutrino and cosmic-ray fluxes due 
to AGNs. Many models for the isotropic neutrino flux from the sum of all 
AGN sources appear in the literature [98] . We consider three models as repre- 
sentative. The flux calculated by Stecker and Salamon [97], labelled AGN-SS 
in Figures 18 and 19, has significant contributions from pp and pry interactions 
in the accretion disk. In the model of Nellen, Mannheim, and Bierman [99], 
labelled AGN-NMB, pp collisions are the dominant neutrino source, leading 
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to a flux 



dN Vl 



dE u 



1.13 xl(T 12 

VI TeV 7 



-2 



cm 



-2 „-l 



sr" 1 GeV" 1 , 



(30) 



for E u ^Ax 10 5 GeV. At higher energies one expects the spectrum to steepen, 
because of the lack of parent protons to produce neutrinos. In our rate esti- 
mates, we use the analytic form (30) up to E v = 10 8 GeV and comment on the 
effect of truncating the neutrino energy spectrum. Szabo and Protheroe [100] 
have extended the model of Stecker and collaborators to include all the impor- 
tant energy-loss mechanisms and computed neutrino production in radio-quiet 
AGNs and in the central regions of radio-loud AGNs. Their model results in 
significantly higher fluxes in the energy range between 1 TeV and 10 3 TeV. 
We take the parametrization 

10- 10 - 5 ( 2 cm- 2 s- 1 sr- 1 GeV" 1 , E v & 10 3 TeV, 

,. < 31 > 

(77 \ — o.o 
I cm" 2 s- 1 sr" 1 GeV" 1 , E v k, 10 3 TeV, 
1 TeV/ 




to represent their hardest spectrum, corresponding to a scaled diffusion coef- 
ficient, b = 1. In the interval 1 TeV ^,E U ^ 10 TeV, this flux is in conflict with 
the upper limit determined by the Frejus Collaboration [101]. This curve is 
labelled AGN-SP. The electron-neutrino fluxes are taken to be one-half of the 
muon-neutrino fluxes. 

All of these fluxes are consistent with the upper limits deduced from horizontal 
air showers by the EAS-TOP Collaboration at Campo Imperatore [102]. For 
10 5 GeV < E v < 10 6 GeV, they infer the "all-flavor" bound (v = v e , v e , z/ M , v^) 

10 6 GeV 

/ ^^<1.5xl0- 8 cm- 2 s- 1 sr- 1 . (32) 
J dE v 

10 5 GeV 



Assuming that the spectrum in this interval is proportional to E 2 , they obtain 
a bound on the differential flux, 

dN v , r __ q / E, 



dE u 



< 1.5 x 10- 9 (pj^v) cm_2 s_1 sr_1 GeV" 1 . (33) 



At the v e e — * W resonance energy, the limit on the v e flux is 

< 7.6 x 10- 18 cm- 2 s- 1 sr- 1 GeV" 1 . (34) 
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The remaining curves represent fluxes from two models of neutrino production 
in interactions of cosmic rays with the microwave background photons. These 
fluxes, calculated numerically by Yoshida and Teshima by Monte Carlo meth- 
ods [103], update earlier analytical results [104,105]. The fluxes depend on the 
redshifts of the cosmic-ray sources: the CR-4 flux corresponds to a maximum 
redshift, or turn-on time, of z max = 4 and evolution parameter m — 4, while 
the CR-2 curve corresponds to £ max = 2 and m — 0. The two models represent 
the extremes presented by Yoshida and Teshima. Separate calculations were 
made for electron and muon neutrinos. 



8.2 v a and v a Interactions 



With these representative fluxes, we turn to the calculation of event rates and 
the implications of the new cross sections presented in §3. As we have noted 
in §1, the effective volume of a detector may be considerably enhanced over 
the instrumented volume by recording charged-current v^N interactions that 
occur in the rock or ice surrounding the detector. The upward muon event 
rate is shielded from the flux of atmospheric muons, and has the advantage 
of utilizing more underground target material. Muons produced with = 
10 TeV will travel, on average, a few kilometers as their energy is degraded 
to 1 TeV. The upward muon event rate depends on the v^N cross section 
in two ways: through the interaction length which governs the attenuation of 
the neutrino flux due to interactions in the Earth, and through the probability 
that the neutrino converts to a muon energetic enough to arrive at the detector 
with larger than the threshold energy E™ n . 

For the case of isotropic fluxes, such as the AGN and cosmic neutrino fluxes 
presented in §8.1, the attenuation can be represented by a shadow factor that 
is equivalent to the effective solid angle for upward muons, divided by 2n: 

S{E V ) = [d cos 9 [ d(j)exp [-z(9)/C- mt (E u )] . (35) 



The interaction length £- mt {E u ) is shown in Figures 11 and 12 for vN and 
vN interactions, respectively. The column depth z(0) is plotted in Figure 15. 
We show the shadow factors computed with the CTEQ-DIS, D_ and CTEQ- 
DLA total cross sections in Figure 20. All of these lead to greater shadowing 
than the EHLQ-DLA distributions used in earlier work. In fact, neither the 
charged-current cross section nor the total cross section is quite appropriate in 
the shadow factor. Neutral-current interactions degrade the neutrino energy, 
but do not remove neutrinos from the beam. A full accounting of the effect 
of neutral currents on the underground upward neutrino flux has been given 
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Fig. 20. The shadow factor S(E U ) for upward-going neutrinos assuming that a = <7 to t 
in (35) for CTEQ-DIS (solid line), CTEQ-DLA (dot-dashed) and D_ (dotted) par- 
ton distribution functions. Also shown is the shadow factor using the EHLQ cross 
sections (dashed line). 

in Ref. [105]. We compute the shadow factor using the interaction lengths for 
the charged-current and charged- current plus neutral-current interactions to 
bracket the number of events for a given model. The longer charged-current 
interaction length leads to higher event rates. 

The probability that a muon produced in a charged-current interaction arrives 
in a detector with an energy above the muon energy threshold E™ m depends 
on the average range (R) of a muon in rock, 

1-Ef n /E v 

(R(E„; < n )> = -4ttt / dyR(E u (l - y), E™) ^jf" V) .(36) 

The range R of an energetic muon follows from the energy-loss relation 

- dEjdx = a{E„) + b{E^. (37) 

If the coefficients a and b are independent of energy, then 

R(E„, Ef n ) = \ In a + b ^ . (38) 
b a + bE™ m v ; 

In our calculations below, we use a = 2.0 x 10~ 3 GeV cmwe -1 and b = 3.9 x 
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Fig. 21. Mean ranges in rock on muons produced in charged-current interactions of 
neutrinos with energy E v . The Lipari-Stanev (solid) and analytic (dashed) ranges 
are shown for E" 1111 = 1 and 10 TeV. 

10~ 6 cmwe -1 in this analytic range formula [106]. We have also considered 
muon ranges evaluated numerically by Lipari and Stanev, which include the 
energy dependence in a and b [107]. In Figure 21, we compare the Lipari- 
Stanev (LS) range and the analytic range for E™ m = 1 and 10 TeV. The 
average range is essentially independent of the parton distribution functions, 
as they all have the same general form for da/dy. 

The probability that a neutrino of energy E v produces an observable muon is 
P,(E U , Ef») = N A a cc (E u )(R(E u ; E™)), (39) 

where Na is Avogadro's number. The event rate for a detector with effective 
area A is 

r dN 
Rate = Aj dE v P„(E u] E™»)S{E V )—. (40) 

The — > n probabilities are plotted in Figure 22 for the three new parton 
distribution functions, as well as the EHLQ-DLA parton distributions, for 
E™ in = 1 and 10 TeV. The effect of the larger cross sections is to increase 
the probability that a neutrino produces an observable muon, but also to in- 
crease the attenuation of neutrinos en route to the detector. The net effect is 
that for the CTEQ-DIS, CTEQ-DLA and D_ cross sections, the combination 
P il (E V) E r ^ in )S(E v ) has little dependence on the choice of parton distribution 
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Fig. 22. Probability that a neutrino of energy E v produces an observable muon with 
energy exceeding E™ m = 1 and 10 TeV, calculated for the Lipari-Stanev range. The 
curves correspond to the CTEQ-DIS (solid), CTEQ-DLA (dot-dashed), D_ (dotted), 
and EHLQ-DLA (dashed) parton distributions. 



0.01 




Fig. 23. The product P^(E U , E I ^ m )S(E l/ ), calculated using the Lipari-Stanev range 
and shadow factor determined by the total cross section, for E™ m = 1 and 10 TeV. 
The curves correspond to the CTEQ-DIS (solid), CTEQ-DLA (dot-dashed), D_ 
(dotted), and EHLQ-DLA (dashed) parton distributions. 



functions, as seen in Figure 23. The CTEQ-DLA and D_ distributions yield 
upward event rates within a few percent of those calculated for CTEQ-DIS dis- 
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Table 4 

Upward fi + + fj,~ event rates per steradian per year arising from u^N and u^N 
interactions in rock, for a detector with effective area A = 0.1 km 2 and muon energy 
threshold £™ m = 1 TeV. The smaller value of each pair corresponds to attenuation 
by the total cross section; the larger to attenuation by charged-current interactions. 



CTEQ-DIS EHLQ-DLA 
analytic LS LS 



ATM [92] 


170- 


-173 


138- 


-141 


124- 


-126 


AGN-SS [97] 


106- 


-126 


77- 


-92 


70- 


-82 


AGN-NMB [99] 


134- 


-146 


102- 


-111 


93- 


100 


AGN-SP [100] 


3570- 


-3870 


2740- 


-2960 


2440- 


-2660 



Table 5 

Upward /i + + fi~ event rates per steradian per year arising from u^N and u^N 
interactions in rock, for a detector with effective area A = 0.1 km 2 and muon 
energy threshold E™ m = 10 TeV. The smaller value of each pair corresponds to 
attenuation by the total cross section; the larger to attenuation by charged-current 
interactions. 



Flux 


CTEQ-DIS 
analytic LS 


EHLQ-DLA 
LS 


ATM [92] 


4 


3 


3 


AGN-SS [97] 


62-75 


43-51 


39-46 


AGN-NMB [99] 


42-49 


30-34 


27-31 


AGN-SP [100] 


1060-1200 


747-843 


683-760 



tributions for E™ m = 1 TeV. Consequently we only show rates corresponding 
to the EHLQ and CTEQ-DIS cross sections. 

In Table 4 we show the upward-muon event rates for a detector with an ef- 
fective area of 0.1 km and a muon energy threshold of 1 TeV. These event 
rates are for muons and antimuons with modern (CTEQ-DIS) and ancient 
(EHLQ-DLA) parton distribution functions, and show the difference between 
the analytic and Lipari-Stanev muon ranges. As a practical matter, we have 
taken the upper limit of the energy integral (40) to be E™ x = 10 8 GeV, the 
limit of the Lipari-Stanev analysis of the muon range. The event rates from 
atmospheric neutrinos are roughly comparable to the AGN neutrino event 
rates for this muon energy threshold. In fact, most of the AGN event rate 
comes from the first few energy decades. The ATM rate comes entirely from 
E v < 10 6 GeV. For the AGN-SS flux, only about 75% of the rate comes from 
E v < 10 6 GeV, but by E v = 10 7 GeV, one has essentially all of the rate. About 
5% of the AGN-NMB rate comes from neutrino energies above 10 6 GeV. The 
details of the turnover of the AGN-NMB energy spectrum should not affect 
the predicted event rate significantly. The rates calculated with the CTEQ- 
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DLA and D_ distributions give numerical rates essentially indistinguishable 
from the CTEQ-DIS numbers, as one would expect from Figure 22. 

Table 5 shows the upward /i + + jjT event rate for a muon energy threshold of 
10 TeV. The atmospheric neutrino background is significantly reduced. About 
85% of the AGN-NMB rate arises from neutrino energies below 10 6 GeV, 
an indication that the expected steepening of the spectrum may reduce the 
event rates reported in the table by some 10 to 20%. In the AGN-SS model, 
for which the spectrum is predicted beyond E v = 10 9 GeV, more than 95% 
of the rate comes from E v < 10 7 GeV. Integrated over 2tt solid angle, the 
annual rates are very encouraging. We expect at least 190 AGN events on 
a ten-percent background. The atmospheric-muon background is negligible, 
except at the Earth's surface, where horizontal muons must be avoided. If the 
Szabo-Protheroe fluxes are correct, contrary to the Frejus evidence [101], the 
detection of diffuse astrophysical neutrinos is imminent. 

The cosmic-neutrino fluxes shown in Figure 18 are of interest for E v > 10 7 GeV. 
To evaluate the event rate for cosmic-neutrino interactions, we have evaluated 
the energy integrals from 10 7 GeV to 10 12 GeV using the analytic formula 
for the muon range. The upward event rates for muons with energies above 
10 7 GeV are shown for a variety of parton distributions and detector condi- 
tions in Table 6. The CTEQ-DIS cross sections yield upward rates only about 
20% larger than those implied by the EHLQ-DLA cross sections. The upward 
muon event rate appears to be very difficult to observe in a 0.1-km 2 detector. 

To further explore the possibility of detecting cosmic neutrinos, we turn our 
attention to the downward and horizontal v^N event rates. The passage of 
neutrinos through the Earth reduces the upward angle-averaged neutrino flux 
by a factor of ten at E v = 10 7 GeV, and even further as the neutrino energy 
increases. The cosmic neutrino energy spectrum is nearly flat for E u between 
10 7 GeV &E V & 10 9 GeV, so that in the absence of shadowing, the rate would 
be dominated by neutrinos with energies near the upper end of that range. 

It is a good approximation to set the shadow factor to unity when considering 
downward neutrinos and for incident angles such that the column depth z 
of the intervening rock is small compared to the neutrino interaction length. 
The range of interaction lengths for E v = 10 7 - 10 9 GeV is 6.6 x 10 3 - 9.4 x 
10 2 kmwe. For a detector at the surface of the Earth, these lengths correspond 
to angles between 1.3° and 8.9° below horizontal. The detectors of interest 
are kilometers underground, so the precise angle at which the column depth 
equals the interaction length depends on details of the location of the detector. 
However, the numbers indicate in general that one can reliably set S(E U ) = 1 
only for neutrinos that are entering the detector from above or horizontally. 

The downward event rates in Table 6 are calculated with no shadowing. Two 
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sets of downward rates are shown: the first is for contained events in an effective 
volume of V e d = A ■ 1 km = 0.1 km 3 for E v between 10 7 and 10 12 GeV, while 
the second set corresponds to V c s = A(R), with E™ m = 10 7 GeV. 

For the contained events, the downward muon event rate is enhanced relative 
to the upward rate by a factor of 3 to 4 for the CTEQ-DIS parton distribution 
functions. Differences in parton distribution functions are much more striking 
in the downward event rate than for the upward event rate. In this case, the 
CTEQ-DIS rates are about twice as large as the old EHLQ-DLA rates. Even 
with the most optimistic flux and the highest (D_) estimate of the neutrino- 
nucleon cross section, the contained rates for cosmic-neutrino interactions are 
very low. 

A larger rate of muons from cosmic neutrinos would obtain if it were feasible to 
take advantage of the average muon range of about 10 km. The second set of 
downward event rates uses the analytic range to establish the effective volume. 
The location of the detector will limit the range enhancement of the effective 
volume, since none of the planned detectors will be deployed at a depth of 
10 km. Even if one could take advantage of the full range enhancement over 
2n solid angle, the predicted rate for the CR-4 flux using the D_ cross section 
is on the order of 0.3 event per year in a detector with A = 0.1 km 2 . 

In our discussion of the downward event rates, we have not addressed the 
problem of the atmospheric muon background. For uncontained events, the 
neutrinos must interact to produce a muon signal, while the muons produced 
in the atmosphere by cosmic rays need only pass through the detector volume 
to be recorded. At = 10 7 GeV, the flux of muons is comparable to the 
flux of neutrinos in the CR-4 model at the Earth's surface [94]. Underground, 
the muon energy is degraded according to the range formula (37). To a good 
approximation, the muon flux is decreased by a factor of exp(— bqz) when 
dN/dE^ oc E~( 1+l \ The vertical muon flux of (29) corresponds to 7 = 2.7. 
Taking b = 3.9 x 10 -6 cmwe -1 as before, we find that the energy spectrum 
of the atmospheric muons below ground is degraded by a factor of about 
exp(— l.lz krnwe" 1 ). At a column depth of 8 kmwe, the suppression amounts 
to a factor of ~ 10~ 4 . Since the neutrino-to-muon conversion rate involves 
the multiplicative factor N A <j cc (E v )L « 10-\E V /10 7 GeV) a4 (L kmwe -1 ), 
the background from atmospheric muons is a concern at depths substantially 
less than 8 kmwe. To compensate, the solid angle must be restricted to include 
only large column depths. Consequently, it is overly optimistic to assume that 
uncontained neutrino-induced events can be observed over a 2n solid angle. 
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Table 6 

The u~ + u + event rates per steradian per year corresponding to two models of 
the cosmic neutrino flux (CR-2 and CR-4 [103]), for a detector with effective area 
A = 0.1 km 2 and muon energy threshold E™ n = 10 7 GeV. For upward events, we 
calculate the attenuation using the total cross section. For downward events we set 
S(E V ) = 1. 



Effective Volume 


direction 


parton distributions 


CR-2 




CR-4 




A 


/ T~>\ 

■ (R) 


upward 


CTEQ-DIS 


1.9 x 


10" 


-5 


1.1 x 10- 


-3 


A 


■(R) 


upward 


L/ 1 r^Q-DLA 


1.8 x 


10~ 


-5 


1 n w i n- 
l.U X 1U 


-3 


A 


■(R) 


upward 


D_ 


1.9 x 


10~ 


-5 


1.1 x 10- 


-3 


A 


■(R) 


upward 


EHLQ-DLA 


1.6 x 


io- 


-5 


9.2 x 10- 


-4 


A ■ 


1 km 


downward 


CTEQ-DIS 


7.4 x 


io- 


-5 


3.5 x 10- 


-3 


A ■ 


1 km 


downward 


CTEQ-DLA 


5.1 x 


io- 


-5 


2.6 x 10- 


-3 


A ■ 


1 km 


downward 


D_ 


1.2 x 


io- 


-4 


4.9 x 10- 


-3 


A ■ 


1 km 


downward 


EHLQ-DLA 


3.4 x 


io- 


-5 


1.8 x 10- 


-3 


A 


■(R) 


downward 


CTEQ-DIS 


1.0 x 


io- 


-3 


3.3 x 10- 


-2 


A 


■(R) 


downward 


CTEQ-DLA 


7.1 x 


io- 


-4 


2.5 x 10- 


-2 


A 


■(R) 


downward 


D_ 


1.7 x 


io- 


-3 


4.8 x 10" 


-2 


A 


■(R) 


downward 


EHLQ-DLA 


4.7 x 


10- 


-4 


1.7 x 10- 


-2 



13 v e and v e Interactions 



Finally we turn to the calculation of event rates involving electron neutrinos. 
Calculations for v e N charged- current interaction event rates proceed as above 
with u^, except that the electron range is significantly shorter than the muon 
range. In general, only contained events can be observed because of the rapid 
energy loss (or annihilation) of electrons and positrons. Accordingly, event 
rates for electron neutrinos are smaller than muon event rates by the flux ra- 
tio times the detector length divided by the mean muon range. However, the 
rapid development of electromagnetic showers may make it possible to detect 
upward-going air showers initiated by an electron neutrino that interacts near 
the surface of the Earth. The Landau-Pomeranchuk-Migdal effect [108,109] 
enhances the distance an electron can travel in the Earth. For electrons pro- 
duced in u e N interactions at energy E u , the mean path length is 



Llpm(E„) ~ 40 cmwe 



(i - (y(K))) 



62 TeV 



1/2 



(41) 
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Table 7 

Downward resonant v e — > VF~ events per steradian per year for a detector with 
effective volume V e g = 1 km 3 . Also shown are the potential downward (upward) 
background rates from v^N and v^N interactions above 3 PeV. 



Flux 


v e e -> P^fx 


v e e — > hadrons 


(*V P„)JV CC 


P M )iV NC 


AGN-SS [97] 


6 


41 


33 (7) 


13 (3) 


AGN-SP [100] 


3 


19 


19 (4) 


7(1) 



A very-large-area air shower array might therefore constitute a large-volume 
detector for electron neutrinos. 

Prospects for the detection of electron antineutrinos are more favorable around 
6.3 PeV, the energy for resonant W~ formation in v e e collisions. The contained 
event rate for resonant W production is 

10 (^rv)V2m 
Rate = _ N A V eS J d E Pe a Pee (Ep e )^. (42) 

{M w -2V w ) 2 /2m 



We show in Table 7 the number of resonant v e e events produced per steradian 
per year in a 1-km 3 detector for two models of the diffuse neutrino flux from 
AGNs that apply in this energy regime. We recall that, at the resonance en- 
ergy, upward-moving electron antineutrinos do not survive passage through the 
Earth. The form oc (1— y) 2 of the differential cross section (20) for v e e — > z^u - 
means that the mean energy of muons arising from W~ formation and decay 
will be (Efj) ~ \E™ S wl.4 PeV. The resonance signal is not background-free. 
We have also gathered in Table 7 the downward and upward rates for the 
charged- current iy^N — > ii~ + anything and v^N — > ji + + anything) back- 
ground to the v e e — > W~ — > signal, and the downward and upward rates 
for the neutral-current [y^N — > + anything and v^N — > P M + anything) 
background to the ^ e e — >■ — > hadrons signal. For this background estimate 
we have included all events induced by neutrinos with energies above 3 PeV. 
At the surface of the Earth, (29) leads to an estimate of 5 atmospheric-muon 
events per steradian per year above 3 PeV. Better discrimination against back- 
ground is clearly desirable. 



9 Summary and Outlook 



We have studied the implications of new knowledge of nucleon structure at 
small values of x for the detection of ultrahigh-energy neutrinos from extrater- 
restrial sources. Using a variety of modern parton distributions, we have calcu- 
lated cross sections for the charged-current reactions, u^N — > iT + anything 
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and v^N — * /i + + anything, that will be used to detect UHE neutrinos. Up 
to energies of about 10 16 eV, parton distributions that entail different behav- 
iors as x — > yield very similar cross sections. The calculated cross sections 
are in good agreement with the charged-current cross section inferred from 
e~p interactions at HERA at an equivalent neutrino energy of 4.7 x 10 13 eV. 
At energies below 10 15 eV, the new cross sections are about 15% larger than 
those calculated by Quigg, Reno, and Walker using the EHLQ structure func- 
tions and the double logarithmic approximation for the approach to x — 0. At 
higher energies, the difference between new and old cross sections increases 
rapidly, reflecting the HERA observation of large parton densities at small x. 
At 10 20 eV, our nominal cross sections, calculated from the CTEQ3 parton 
distributions, are about 2.4 times the EHLQ-DLA cross sections of a decade 
ago. In the regime above 10 16 eV, the cross sections are sensitive to parton 
distributions at very small values of x, where there are no direct experimental 
constraints. Accordingly, different assumptions about the x — > behavior lead 
to different cross sections. At 10 20 eV, the resulting uncertainty reaches a fac- 
tor of 2 ±x . We have also calculated the neutral-current v^N — > + anything 
and v^N — > P M + anything cross sections that contribute to the attenuation 
of UHE neutrinos as they traverse the Earth. 

We have estimated event rates in large-volume detectors for downward- and 
upward-moving muons produced in charged-current interactions. The increased 
charged-current cross section translates directly into increased downward event 
rates, but the observation of downward events is complicated by the back- 
ground of cosmic-ray muons. For upward events, the increased interaction 
rate is nearly compensated by the increased attenuation of UHE neutrinos in 
the Earth. 

We expect that the new generation of neutrino telescopes will detect UHE neu- 
trinos from extraterrestrial sources, and will begin to test models for neutrino 
production in active galactic nuclei. For the CTEQ-DIS cross sections and the 
Lipari-Stanev muon range-energy relation, we find that in one steradian-year, 
a detector with an active range of 0.02 km 2 would record between 16 and 592 
upward-moving muons with energies above 1 TeV produced by interactions 
of AGN neutrinos, on a background of about 28 events produced by atmo- 
spheric neutrinos. If the muon energy threshold is raised to 10 TeV, the rates 
induced by diffuse AGN neutrinos will be between 9 and 170 events on a back- 
ground of less than one event. The range of signal events reflects the spread 
in predictions of the diffuse neutrino flux from AGNs. 

The outlook for the detection of cosmic neutrinos at energies around 10 17 eV 
is less encouraging. Even in a detector with an effective volume of 1 km 3 , the 
most favorable model for the cosmic-neutrino flux leads to less than one event 
per steradian-year with > 10 16 eV. 



41 



Finally, we have considered the reaction v e e — > W~ as a means of probing the 
v e spectrum in the neighborhood of the resonant energy, 6.3 x 10 15 eV. We 
estimate that a detector with effective volume 0.2 km 3 would record between 
4 and 7 downward v e e — > W~ — > P^/x events and between 24 and 50 down- 
ward hadronic events per year. The backgrounds from deeply inelastic u^N 
scattering are not negligible. 

We are optimistic that progress toward large-volume neutrino telescopes, ini- 
tially based on water- Cerenkov and ice-Cerenkov techniques, will soon lead 
to the detection of ultrahigh-energy neutrinos from extraterrestrial sources. 
With the ability to detect UHE neutrinos will come the possibility of looking 
deep within some of the most energetic structures in the universe. For neutrino 
energies up to 10 16 eV, which spans the range of interest for testing models 
of active galactic nuclei, the neutrino-nucleon cross sections can be predicted 
with confidence. We expect neutrino telescopes to emerge as an important 
astrophysical tool. 
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